LATTICE STABILITY FOR ATOMISTIC CHAINS MODELED BY LOCAL 
APPROXIMATIONS OF THE EMBEDDED ATOM METHOD 



XINGJIE HELEN LI AND MITCHELL LUSKIN 

Abstract. The accurate approximation of critical strains for lattice instability is a key criterion 
for predictive computational modeling of materials. In this paper, we present a comparison of the 
lattice stability for atomistic chains modeled by the embedded atom method (EAM) with their 
approximation by local Cauchy-Born models. We find that both the volume-based local model 
and the reconstruction-based local model can give 0(1) errors for the critical strain since the 
embedding energy density is generally strictly convex. The critical strain predicted by the volume- 
based model is always larger than that predicted by the atomistic model, but the critical strain for 
reconstruction-based models can be either larger or smaller than that predicted by the atomistic 
model. 



1. Introduction 

Predictive multiscale computational methods must be accurate near lattice instabilities that 
characterize the formation and movement of cracks, dislocations, and grain boundaries. In this 
paper, we present analytic results comparing the lattice instabilities predicted by an atomistic 
chain modeled by the embedded atom method (EAM) with the lattice instabilities predicted by 
local approximations of the atomistic model. 

Since it is not possible to compute large enough fully atomistic systems to accurately approxi- 
mate the interaction of local defects with long-range elastic fields, atomistic-to-continuum coupling 
methods have been proposed [H[^[51 [TT1 - [I3t [TB iri8ll2Uj . For crystalline solids, the continuum region 
is generally computed by coarse-graining a local approximation of the nonlocal atomistic model. 
An atom in the nonlocal atomistic region interacts with all of its neighbors within a cutoff radius. 
In the continuum region, the Cauchy-Born rule is used to derive a local model that approximates 
the interactions of atoms beyond their nearest neighbors by modified interactions of their nearest 
neighbors. 

To verify that an atomistic-to-continuum coupling method accurately reproduces the lattice 
stability of the fully atomistic model, it is necessary to first verify that the local (continuum) model 
itself reproduces the lattice stability of the fully atomistic model. Even if the local model reproduces 
the lattice stability of the fully atomistic model, the atomistic-to-continuum coupling method may 
not reproduce the lattice stability of the fully atomistic model because of the error introduced by 
the coupling [3|[T9]. 

It has been proven in [3] that the Cauchy-Born local model reproduces the lattice stability for an 
atomistic chain modeled by Lennard-Jones type pair interaction (we note that the volume-based 
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and reconstruction-based local models are equal for a pair potential interaction). And it has been 
proven for multidimensional lattices that the set of stable uniform strains for the atomistic model 
is a subset of the set of uniform strains for the Cauchy-Born volume-based local model [ElE], but 
the equality of these sets has not been demonstrated analytically. As a matter of fact, numerical 
experiments in [8] suggest that the the inclusion is strict in some cases. 

In this paper, we prove for an atomistic chain that not only are the sets of stable uniform strains 
different for the atomistic model and the local models for a many-body potential, but the set of 
stable uniform strains can be different for volume-based and reconstruction-based local models. We 
will focus our analysis on the embedded atom method [ZllHKUj, which is an empirical many-body 
potential that is widely used to model FCC metals such as copper and aluminum. We identify the 
critical assumptions for the pair potential, electron density function, and embedding function to 
study the lattice stability of the atomistic and the different local models. We find that both the 
volume-based local model and the reconstruction-based local model can give 0(1) errors for the 
critical strain since the embedding energy density is generally strictly convex. 

In Section [21 we present the notation used in this paper. We define the displacement space U 
and the deformation space yp. We then introduce the norms we will use to estimate the modeling 
error and the displacement gradient error. In Section [3l we briefly review the formulae of the fully 
atomistic EAM model and the volume-based and the reconstruction-based local quasicontinuum 
(QCL) model, respectively. 

In Section m we give precise stability estimates for the fully atomistic model, the volume-based 
and the reconstruction-based local models for a uniformly strained chain. We then compare the 
stability conditions of each model under different assumptions. We summarize our results and 
discuss extensions to multidimensional issues in the Conclusion. 



2. Notation 

In this section, we present the notation used in this paper. We deflne the scaled reference lattice 

eZ := {el:ie Z}, 

where e > scales the reference atomic spacing and Z is the set of integers. We then deform the 
reference lattice eZ uniformly into the lattice 

FeZ := {Fei : £ G Z} 

where F > is the macroscopic deformation gradient, and we define the corresponding deformation 
yp by 

{yF)i '■= Fel for — oo < ^ < oo. 
For simplicity, we consider the space U of 2A^-periodic zero mean displacements u = (u^)fgz from 
yp given by 

hi := Ivl : ui^2N = Ui for £ £ Z, and = >, 

i=-N+l ^ 

and we thus admit deformations y from the space 

yp := {y : y = yp + u for some u £ U}. 

We set € = 1/N throughout so that the reference length of the periodic domain is fixed. 
We define the discrete differentiation operator, Du, on periodic displacements by 

(iJujf := , — oo < £ < oo. 
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We note that (Du)^ is also 2iV-periodic in i and satisfies the zero mean condition. We will denote 
{Du)^ by Du£. We then define 

y)u)^:=^^^^^f^, -oo<£<oo, 

and we define {D^'^^u)^ and {D^'^^u)^ in a similar way. To make the formulas concise and more 
readable, we sometimes denote Dui by u'^, D^'^^ui by u'^, etc., when there is no confusion in the 
expressions. 

For a displacement u and its discrete derivatives, we define the discrete £^ norms by 

/TV \ / ^ \^^'^ 

l|u||^2:=le ^ > l|u'||^2:=le ^ Kl"^ \ , etc. 

\ e=-.N+i / \ e=~N+i / 

Finally, for smooth real- valued functions £{y) defined for y € yp, we define the first and second 
derivatives (variations) by 



{5£{y),w) := T:—{y)we for ah G ZY, 



l=-N+l 

{5'^£{y)v,w) := ^ — — — {y)viWm for ah v, w G Z^. 

i,m=~N+i "^y^^y^ 

3. The Embedded Atom Model and Its Local Approximations. 

In this section, we will give a short description for the next-nearest neighbor atomistic EAM 
model and its approximations . 

3.1. The Next-Nearest Neighbor Atomistic EAM Model. Given deformations y G yp, the 
total energy per period of the next-nearest neighbor atomistic EAM model is 

Ci(y):=^^"(y)+-^(y), (3.i) 

where £"'{y) is the total atomistic energy and J-{y) is the total external potential energy. The total 
atomistic energy £°'{y) is the sum of the embedding energy, £°'{y), and the pair potential energy, 
£"'{y). The energy expression is 

N 

£\y) := t{y) + £\y) = e (^/(y) + ^"(y)) • (3-2) 

^=-Ar+l 

The embedding energy per atom (per atomistic reference spacing e) is defined as £1{y) ■= G (yo"(y)) , 
where G{p) is the embedding energy function and p'}{y) is the total electron density at atom £: 

Peiy) ■= p{y'i) + p{v\ + + p(%+i) + + %+2)- 

The function p{r/e) is the electron density contributed by an atom at distance r. 
The pair potential energy per atom (per atomistic reference spacing e) is 
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where cj){r/e) is the pair potential interaction energy [7]. Our formulation allows general nonlinear 
external potential energies J-{y) defined for y G yp, but for simplicity, we only consider the total 
external potential energy for 2A^-periodic dead loads f 

N 

^{y) ■=- X] ^f^y^- 

i=-N+l 

The equilibrium solution y" of the EAM-atomistic model ()3.ip then satisfies 

-{5£-{y-)^^) = -(5^«(y"),w) - (<5f~'^(y'^),w) = (<5^(y"),w) for ah w G ZY. (3.3) 
Here the negative of the embedding force of (j3.3|) is 

N 

(5^'^(y"),w)=e G'[pny''))-[p'{Dy'l)w', + p'{Dy'l + DyU){< + <_,) 

e.=-N+i 

+ p'iDy^^^)w'^^, + p'{Dy^^^ + Dy^^^)^^! + ^£+2) 
the negative of the pair potential force of (j3.3|) is given by 

N 

2 



^ 1 r 



e=-N+i 



and the negative of the external force is formulated as 

((5J"(y),w)= —{y)we = - ^ efeWi. 

e=-N+i e=-N+i 

3.2. The Local EAM Approximations. In this subsection, we will briefly review the idea of 
the two different local approximations, the volume-based and the reconstruction-based, and give 
their expressions respectively. 

3.2.1. The Volume-Based Local EAM Approximation. The idea of the volume-based local approxi- 
mation based on the Cauchy-Born rule was first proposed in |12 y i5 ^ [T7]. We denote this energy by 
£'^''"(y), and we can formulate the local energy associated with each atom as 

sriy) ■■= sriy) + sriy) =\g [pTm + \g (Pl^i(y)) 

+ \ We) + mv'i) + <\>h)'i+x) + '^(2?/;+i)] , 
where the total local electron density at atom i is 

(y) := 2p(y;) + 2p(2y;). 
Then the total volume-based local energy is 

AT TV 

The equilibrium solution y'^'" then satisfies 

-(5£:'=''^(y^'^), w) = -(^^'^.-(y-'''), w) - l,bt^\y-^^),^) = (5^(y^'^'), w) for ahw G U. (3.5) 
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The negative of the embedding force of ()3.5p is 

N 

{5£^'^{y^n,^) =e i^' (/^^(y''")) " [p'iDyr) + 2p'{2Dyr)] w', 

£=-Af+l 

and the negative of the pair potential force of p.Sp is given by 

^ 1 

e=-N+i 

3.2.2. The Reconstruction-Based Local EAM Approximation. Using the Cauchy-Born approxima- 
tion, one can also reconstruct the position of each atom [S] and compute the energy £^''^{y) by the 
approximation 

sriy) = s^iy) + sriy) =g (p'/iy)) + \ Wi) + H2y'i) + M+i) + mWi)] > 

where the reconstruction-based local electron density at atom I is 

pYiy) ■= piy'e) + Pi^y'e) + Pive+i) + P(2y^+i). 
Thus, the total energy of the reconstruction-based local model is 

Stoliy) {y)+£'''{y) + Hy) 

N 

=e ^ [G[p{y',) + p{2y',) + p{y',^,) + p{2y',^,)] 

e-N+i (3-6) 

1 ^ 

+ 2 Wi) + ^C^y'i) + ^iy'i+i) + ^i2y',^i)] } - e E f^y^- 

e=-N+i 

The volume-based and reconstruction-based local energies have the same pair potential energy, but 
their approximations for the embedding energy are quite different. 

We compute the equilibrium solution of the reconstruction-based local model ()3.6p from 

-{SS^'^iy^n,^) = -{Si'''{y'n,^) - {S£'''{y'n,^) = {SHy'n,^) for ah w G ZY. (3.7) 
Here the negative of the embedding force of (|3.7p is 

{5t'^{y^n,^) =e J2 (p'/iy'^n) ■ [{p'iDyr) + 2fl{2DyY)) w[ 

l=-N+l 

+ {p'iDyl'l^) + 2p'{2Dyl'l^)) w[^^] , 
and the negative of the pair potential force of (j3.7p is 

^ 1 

{5£^^^{y^n,^) =e E 2 ^ WiDyY) + 2</<'(2Z)yf )] w', + [0'(Z)y,%\) + 2</>'(2Z)y,%^,)] w',^,] . 
e=-N+i 

The pair potential energy of both local approximations are exactly the same, but the embedding 
parts are quite different, which leads to different critical strains for lattice instability. We will 
analyze the lattice stability for all of the models in the next section. 
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4. Sharp Stability Analysis of The Atomistic and Local EAM Models. 

In this section, we analyze and compare the conditions for lattice stability of the atomistic 
model and the two local approximations for the next-nearest neighbor case. We will use techniques 
similar to those presented in [3j for the atomistic and quasicontinuum methods with pair potential 
interaction. 

4.1. Stability of the Atomistic EAM Model. We first consider the fully atomistic model. The 
uniform deformation y_p is an equilibrium of the atomistic model (|3.2p without external force. We 
call yp stable in the atomistic model if and only if {S'^S^'iyp) is positive definite, that is, 

{6^£''{yF)u, u) = {6^£''{yF)u, u) + {6^£''{yF)u, u) > for aU u G \ {0}. (4.1) 

We computed {6'^£"'{yF)u,u) in [3| to obtain 

{6^£''{yF)u,u) = AfWDuWI - eV2Fp^'^u|||, (4.2) 

where 

Af := ^F + for ^'j, := 0"(F) and c^'^f ■= 0"(2i^) (4.3) 

is the continuum elastic modulus for the pair interaction potential. Thus, we focus on (5^i5"(yi?)u, u), 
which can be formulated as 

{5^£''{yF)u, u) = e ^ (Gp [ppiue + u'e+i) + p2Fiue-i + Ui + u'^+i + ^+2)] ^ 

£=-N+l { 

+G'f [p'Mf + pU< + uU? + p"f{uWi? (4.4) 

where we use the simplified notation 

p'f:=p'{F), pI:=p"{F), p'^p := p{2F), p'^^ := p" {2F) , 
G'f := G'ipnyF)) = G'ipriyp)) = Cif/iyF)), 
G'^ := G"{-pt{yF)) = G"{-pr{yF)) = G"{-pY{yF)). 
We define the continuum elastic modulus for the embedding energy to be 

Af := 4G^ {p'f + 2p'^p) ' + 2G'f {p"f + ^p'ip) , (4-5) 

and we define 

-^4z.. Er^ ■= -\d^'^„^G'L((nL)^ ^20(n'r.r.)^ ^^2n'„n'r.r^ ^GL(2n'^„) . 

(4.6) 



Af := Af + Af, Bp := - [^'2'^ + G'^ [{p'f? + 20(p'2p)2 + Up'pp'^pj + G'p {2p'^f) 

Gf := {8{p'2F? + 'ip'Fp2F) , and Dp := -G'^ (^'2^.)' 
Then (14. ip becomes 

{5^£''iyF)n,u) =A^||Z)u||| + e2SF||Z5(2)^||| ^ ^4^^||^{3)^||| ^ ^6^^||^(4)^||| ^ ^4_7) 

where the detailed calculation can be found in the paper |10j . 

We will analyze the stability of {6'^£"'(yp)u,u) by using the Fourier representation [8] 

k=-N+l^ ^ ^ 



(4.8) 
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We exclude A; = since Du must satisfy the mean zero condition ^^_^_|_^ Dui = 0. 
It then follows from the discrete orthogonality of the Fourier basis that 

N ( 



\Ck\ 



Af + Bf 



k=-N+l 



+ Cf 



4 sin" 



4sin"^ 



kTT 

2N 



kn 
2N 



+ Df 



(4.9) 



4 sin'' 



2N 



We see from (j4.9p that the eigenvalues for k = 1, . . . , of (5^i5"(yF)u, u) with respect to the 
1 1 -Dull £2 norm are given by 

\l = \Usk) for Sk=Asin^(^} 



2N J 
where 

X%{s) := Af + Bfs + Cfs^ + Dfs^. 

The energy and electron densities figures in (7j and jl4j satisfy the following conditions which 
we shall assume in our analysis 

> 0, < 0; p'f < 0, p'2F < 0; > 0, p'^f > 0; and > 0. 

We can derive from the assumption (|4.10p that 

Cf > 0, Df < 0, and 8\Df\ < Cf- 

Since (14. IIP implies that \Dfs\ < 4\Df\ < Cf/2, for < s < 4, we have that 



X1,'(s) = Bf + 2Cfs + 3Dfs'^ > Bf + for aU < s < 4. 

We note from (|4.12|) that the condition Bf > or equivalently 

{p'f^ + 20 (p'2f)' + WfP2f\ + G'f 2p'^F = -Bf < 0, 



'■'2F + 



implies that A^(s) is increasing for < s < 4. We thus conclude that if Bf > 0, then 



((5^^:"(yF)u,u) > A^(si)||L'u|||2 > I^Af + AfJ ||-Du||| for aU u G ZY. 
This result is summarized in the following theorem: 



(4.10) 
(4.11) 

(4.12) 

(4.13) 

(4.14) 



Theorem 4.1. Suppose that the hypotheses (j4.10p and Bf > hold. Then the uniform deformation 
Yf is stable for the atomistic model if and only if 



XUs,) = Af + Bf [iW (^)] + Cf [Asm' (^) 



\2N. 

Ap + Af - 4sin2 + 



+ Df 



4sin2f^) 



[p'f^ + 20 {p'^fY + WfP2f] + g'f 2p'^f} 



+ 4^sin^(^) G'^ [v{p'2Fy + 2p'Fp'2F 



We note that the differences between Sk and Sk^i and between A^(sfc) and A^(sfc-i) are of 

order O ^^^r^ = O (2/ce^) for k = 1, . . . ,N. When the number of atoms is sufficiently large, 

mino<s<4 A^(s) can be used to approximate mini<fc<jv A^(sfc) with 1 < k < N with error at most 
of order 0(e) since A^e = 1. 
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When Bp < and is sufficiently large, the minimum eigenvalue of S'^S'^iyp) is no longer 
Xp{si) and is given by the following theorem. 

Theorem 4.2. Suppose that the hypotheses ()4.10p and Bp < hold, and the number of atoms N 
is sufficiently large. Then Xp{si) defined in Theorem \4-l\ will no longer be the minimum eigenvalue 
of the second variation {6'^ S"" {y p)u, u) . Instead, the minimum eigenvalue will be given by Xp{sk*) 
for some s^*, 1 < k* < N , that is either equal to 4 or close to 




-3BpDp 



with difference of order O (2A;*e^) . 

Proof. Here we will briefly discuss the role of the coefficient Bp and leave the rigorous discussion 
of mino<s<4 Xp{s) under the condition Bp < to section O 

The assumption Bp > guarantees that u'^ = sin(e£7r) is the eigenfunction corresponding to the 
minimum eigenvalue of 6'^£"'{yp) with respect to the norm ||L'u||£2. In fact, when Bp < 0, we have 
A^'(O) < and thus X%{0) will be strictly larger than A^(s*). 

We note that the condition Bp > cannot be expected to generally hold for EAM models when 
the nearest neighbor term Gp{p'p)'^ > dominates. We note, however, that generally G'p < for 
F < 1 [H], in which case Bp >0 is more likely to hold for compressive strains F < 1. □ 

Remark 4.1. We would like to point out that when is small, Xp{si) may be still the minimum 
eigenvalue of 5'^£°'{yp) even if Bp < 0. This is because Xp{sk) is defined on the discrete domain 
1 < A; < A^, so the continuous function Xp{s) is not a good approximation unless is sufficiently 
large. 

4.2. Stability of the Volume-Based and the Reconstruction-Based Local EAM Models. 

In this subsection, we will give stability estimations for the volume-based and the reconstruction- 
based local models, respectively. 



4.2.1. Stability of the Volume-Based Local EAM Model. We focus on the stability of the volume- 
based local model under a uniform deformation yp. Using the equilibrium equation ()3.5p . we obtain 
the second variation 6'^£'^''"{yp) for any u G W \ {0} 

{6^ {yp)u,u) = (Ap + Ap) \\Du\\% = Ap\\Du\\l 



where Ap and Ap are defined in (j4.5p and (|4.3p . respectively. It follows that y_F is stable in the 
volume-based local model if and only if Ap := Ap + Ap > 0. We summarize this result in the 
following theorem. 

Theorem 4.3. Suppose that the hypotheses (I4.10p holds. Then the uniform deformation yp is 
stable in the volume-based local model (j3.4p if and only if Ap := Ap + Ap > 0. 



Remark 4.2. Comparing the conclusions in Theorem 14. 1 1 and Theorem 14. 31 we observe that when the 
hypothesis (j4.13p is satisfied, the difference between the minimum eigenvalues of the fully atomistic 
and the volume-based local models is of order O(e^). This result is the same as for the pair potential 
case However, when the assumptions fails, the volume-based local model will be strictly more 
stable than the fully atomistic model, which will be discussed in the next remark. 
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Remark 4.3. The assumption (j4.13p is necessary for the vahdity of Theorem 14.11 We now give an 
exphcit example showing that the uniform deformation can be strictly more stable for the volume- 
based local model (j3.4p than for the fully atomistic model when (I4.13P fails. We consider the 
case 



Then (|4.13p does not hold since it follows from (j4.10p that 

^2F + G"p \[flpf + 20 (p'2f)' + 12p'f/02f] + G'p2(^!,p 

i>2F + G"f [p'f + 2p2f)^ + G'f2p2f + SGp [2 {p2fY + p'fp'2F 



(4.15) 



> 0. 

We define an oscillatory displacement u, corresponding to the k = N eigenmode in the Fourier 
expansion (j4.8p . by 



ue := {-lYe/{2V2). 

Therefore, 

u'e = {-lY/{V2), ||Du||,. = l, u'; = {-lYiV2)/e. 
From (14. 2|) and (14. 4p we can get 



{S'riyF)u,u) = {6'£''{yF)u,u) + {6'£''{yF)u,u) 
^ 1 

= 6 Yl G'^2/,^^ + (c^'; + 40'2V)ll^u|l| + (-e''/'2F)l|i?^'^U 

e=-N+i 

= G'F2f^'F + {41'f + 4</.'2V) - 40'2V = <^F + G'f2(J'f. 



.~,l|2 



(4.16) 



Thus, we can obtain 



inf (52f'^(yF)u, u) < ^'f + G'f2(J'f. 

uGW\{0}, ||Du||(,2=1 



On the other hand, from Theorem 14.31 we have that 



inf (52£:^'''(y^)u, u) = i^+^F = 4 ^'2'^ + G"f {p'f + 2p'2f) + G'^V2F +4>f+G'f2p'^. 

ugw\{o}, ||i:>u||^2=i L ' ' J 



Therefore, from ()4.15p and ()4.16p we have 



inf 

ueU\{0}, ||Du||^2=l 



(J2£:^'^(yi.)u,u) >0^ + Gi.2/,^> 



inf 

ueW\{0}, ||Du|| 2=1 



{6^£''{yF)u,u). 



This inequality indicates that when the assumption (|4.13p fails, the uniform deformation yp can 
be unstable for the atomistic model, but still stable for the volume-based local model. 

4.2.2. Stability of the Reconstruction-Based Local EAM Model. In this case, we do a similar calcu- 
lation for the reconstruction-based local model and derive the second variation 5'^£^'^'{y) from the 
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equilibrium equation given by ([37 

N 

(5V'^(yi.)u,u)=6 {G'HPF + ^P2F)'{< + <-,lf + G'APF + ¥2F) [K)'+K+i)']} 



i=-N+l 



N 



-AT + l 



= AG"p {p'p + 2p'^pf + 2G'p {p'^ + Ap'ip) + ct>'^ + 40'2V 

-e'G"^{p'p + 2p',^f\\D^')vi\\% 
=Af||-Du||| +e2s^||Z)(2)u|||, 

where Ap is defined in (I4.3P and Bp is defined to be 

Bp ■.= -G"p [p'p + 2p'^pf 



I-Dul 



(4.17) 



We recall that for the EAM-atomistic model, the coefficient Bp of ||-Du||?2 in (14. 7p is defined as 



•2F + ({p'f? + 20 + Up'pp'^p] + G'p {2p'ip) 



Bp 



Comparing Bp with Bp defined in ()4.17p . we find that 



'ip + G"p (16 {p'^pf + ^p'pp'^p) + G'^ (2/,^V) 



Bp = Bp - 

The assumption (j4.10p that cj)2p < implies Bp can be positive while Bp is always negative. 
We similarly use the Fourier representation 



Due 



N 

E 

k=-N+l 



V2 



■ exp ( i k—TT 



to analyze the stability of 6'^ £'^''^ {y p) . Again, we exclude k = because of the mean zero condition 
of Du. Prom the discrete orthogonality of the Fourier basis we have 



TV 



{6'£'''^{yp)u,u)= \ck\'-{Ap + Bp 



k=-N+l 
k^O 



(4.18) 



The eigenvalues A^' of {6 £'^'^{yp)u,'a) with respect to the ||-Du||£2 norm are given by 



^k'' = ^"/(sk) iork = l,...,N, 



where 



Sk = 4 sin 



2N 



and Xp^ (s) := Ap + Bps. 



The assumption (j4.10|) implies that Bp <0 always holds, so A^^(,s) is decreasing for < s < 4 and 
the minimum eigenvalue of 6'^£'^'^{yp) is achieved at k = N, i.e, sat = 4: 

min X'^/isk) = X'/{A) = Ap + 4Bp = 2G'p {p"p + Ap'ip) + <l}"p + ^p. 
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The minimum eigenmode is given by the oscihatory displacement = —u'^^j^^ since 

(52£:-.'-(y^)u,u) = [2G'p {p"p + Ap'ip) + + ^'ip] ||Z)u|||. 

We thus have the fohowing stabihty result for the reconstruction-based local model. 

Theorem 4.4. Suppose that the hypotheses (j4.10p holds. Then the uniform deformation yp is 
stable in the reconstruction-based local model (13. 6|) if and only if 

2G'p (p^ + Ap'ip) + (/>^ + A4>'ip > 0. 

Remark 4.4. Comparing the conclusions in Theorem 14.11 Theorem 14.31 and Theorem 14.41 we note 
that when the assumption (j4.13p holds, i.e., Bp > 0, the fully atomistic model is strictly more stable 
than the reconstruction-based local model. The difference between their minimum eigenvalues is 
0(1), not O(e^) as for the volume-based approximation. When the assumption (|4.13p fails, i.e. 
Bp < 0, the conclusion will be different, and we will rigorously analyze this case in section [5j 

5. Comparison of the Stability of the Atomistic and Local EAM Models 

In this section, we would like to give a full discussion of the sharp stability estimates for all of 
the EAM models. Recall that the eigenvalue function of 6^£"'{yF) is 

XUsk) ■■= Ap + Bpsk + Cpsl + Dpsl for s;,. = 4 sin^ ( |^ ) , /c = l,...,iV, 



\2N^ 

where the coefficients Ap, Bp, Cp and Dp are given in the equation (14. 6p . 

To simplify the following analyses, the number of atoms is assumed to be sufficiently large. 
Thus, we use the global minimum of the continuous function Xp{s) := Ap + Bps + Cps"^ + Dps^ 
for < s < 4 to approximate mini<fc<7v A^(sfc). We note that their difference is at most of order 
O {2ke^) < 0(e). 

We recall that 

min A^(s) = A|(0) if Bp > 0. (5.1) 

0<s<4 

To find mino<s<4 A^(s) when Bp < 0, we first evaluate Xpis) at s = 0,4: 

A^(0) =Ap = 4G'; {p'p + 2p'2p)' + 2G'p {p'^ + 4p'^p) + ^'p + icP'^p, 
A^(4) =(I)"p + 2G'pp"f. 
We next compute the first and second derivatives of Xp{s), which are 

A^'(s) =Bp + 2Cps + 3Dps^, 
A^"(s) =2Cp + 6Dps. 

Since Xp'{s) is a quadratic function, we thus have two critical points of Xp{s) when the coefficients 
satisfy 

C| - WpDp > or equivalents (t)'^p + 2G'pp'^p < -G'^ {p'p - 2p'2pf . 

3 

We can summarize the case when Bp <0 and Gj, - SBpDp < by 

min X%{s) = A^(4) if Bp <0 and C| - SBpDp < 0. (5.2) 

0<s<4 



In the case Gp — 3BpDp > , the critical points are 



Cp - JCl - SBpDp Gp + JGj, - WpDp 

Si = — and So = . 

^ -3Dp ^ -3Dp 
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Since Dp < 0, Xp{s) will then have a local minimum at s* = si and a local maximum at S2- The 
corresponding local minimum value is 

A^(s*) =Af + Bps* + CF{s*f + DF{s*f 

=Af + s [-y- — (') 

where we use Xp'{s*) = to get the last equality. We can thus summarize all of the cases by 
min XUs) = Xl(0) if Bp > 0, 

0<s<4 

^min^A^(s) = A^(4) ii Bp < and C], - 3BpDp < 0, (5 3) 

min XUs) = mm{Xl(s*), XIU)} iiBp <0 and Cl - 3BpDp > 0. 

0<s<4 

We note that the minimum eigenvalues of the volume-based and the reconstruction-based local 
models are separately given by the following expressions 

min X'^is) =A^(0) = Ap = AG'^ {p'p + 2p'^pf + 2G'p (p^ + 4p^V) + + ^<t^2F. 
min A^^s) =A=/(4) = Ap - {p'p + 2p'^pf = 2G'p {p"p + Ap'^p) + (p'^ + Act>'ip. 

0<s<4 

5.1. The Volume-based Local EAM Model versus the Fully Atomistic EAM Model. We 

first compare the minimum eigenvalues of the volume-based local and the fully atomistic models. 
Combining the results of Theorem 14.11 and Theorem 14.2^ we have the following theorem. 

Theorem 5.1. The relation of the stability of the volume-based local model and the fully atomistic 
model depends on the sign of 



Bp:-- 



2F + G"p {p'pf + 2Q{p'^pf + I2p'pp'^p + G'p [2p'ip) 



and can be summarized as follows: 

min X%(s) =A^(0) = A^/ if Bp > 0, 

min Xlis) =min{A^(s*),A^(4)} < A^(0) = min A''/ if Bp < 0. 

0<s<4 L y r\ /J r\ ^ 0<s<4 ^ 

This observation indicates that the set of stable uniform strains for the volume-based local model 
always includes that for the fully atomistic EAM model. 

5.2. The Reconstruction-based Local EAM Model versus the Fully Atomistic EAM 

Model. The relation of the minimum eigenvalues for 6^£^{yp) and E^'''' [y p) is more complicated. 
We note that assumption (j4.10p implies 

A^(0) =Ap = AG"p {p'p + 2p'^pf + 2G'p {fl'p + ApU^p) + 0^ + 4</>^V 
>2G'p {(J'p + 4p'2V) + -/-F + = mm^ A^"(s), 



and we have 



A^(4) - mm^ A^^ (s) = -4 {c^!,p + G'p ■ 2p'^p) . 



We thus conclude that if 

^'^p + G'p ■ 2p'ip < 0, (5.5) 
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then 

A^(4) > mill A"/(s). 

The equal sign is achieved if and only if (j)'2p + G'p ■ ^p'^p = 0. We also have the identity 

ct^'ip + G'f ■ 1P2F = -Bf - [g"f {{p'f? + 20(p2f)' + 12pf/02f) • 
We next compare A^(s*) and mino<s<4 Ap'^(s). The difference of these two is 
X%{s*)- min A'/(s) 

0<s<4 

= AG"f{p'f + 2p',fY + s* [^-^(sr) 
= 4G';(p^ + 2p'2^)' 

Cf - ^Jcl - WfDf 6BfDf - 2Cl + Cf (Cf + - 3BfDf 
-3Df 9D^ 

F KHF -r H2FJ ^ _g^^ ^ 27L>2 



F 



-9Df' 

According to the assumption (|5.5p . we can use the definition of Ci? and Df (|4.6p to get 

Bf = - [(P'^f + G'f ■ 2p'iF + G'^ {{p'f? + 12p'fP2F + ^Qip'^F?)] 
> - G"f {{p'f? + 12p'fp'2f + 2Q{p'2F?) . 
We thus can obtain from the above inequality and the assumption ()4.10p that 
\%{s*)- min A^''(s) 



0<s<4 



>4GMp'f + 2p',^)V^ 

> ^G"f {p'f + 2p'2Ff + ^^^'^ ^ ^^^^ ^ ^^^^^^^'^ ^ ^^{P2F?) 

~^P2F 

Therefore, we have that 

min XUs) = min{A^(0), A^(s*), A^(4)} = min X''/{s) if ^'^f + G'f ■ 2p'^F = 0, 

0<s<4 0<s<4 

min X^s) = min{A^(0), A^(s*), A^(4)} > min X'jf{s) if <))'^f + G'f ■ 2p'ip < 0. 

0<s<4 0<s<4 

Now let us turn to the case that the assumption (|5.5p fails, which means 

<l)'iF + g'f ■ 2f)2F > 0. 
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In this case we have the opposite conclusion that the fully atomistic model £°'{y) is strictly less 
stable than the reconstruction-based local model £^''^{y). From the condition (jS.Sp . we have 

A^(4) - min A^^s) = - 4 {c/^'^p + G'p ■ 2p'^p) < 0, 

0<s<4 

i.e., A^(4) < min A'/(s). 
Before comparing Xp{s*) and mino<s<4 A^^(s), we recall that s* exists if and only if 

(l)2F + 2G'pP2F <^Gp [p'p - 2p'2pf . 

Thus, we actually consider the case that 

< (b'^p + 2G'pp'^p < ^G'^ {p'p - 2p'2Ff . 

We substitute Bp = - [4>2p + G'p ■ 2p'^p + Gp {{p'pf + l2p'pp'2F + 20(^2^)^)] into the inequality 
()5.6p and get 

A^(s*)- min A"/(s) 

>AG"p{p'p + 2p'2pf + ^^ 

Ar" M ^9 ' ^2 , l <P'iF + G'F-2p'ip , G"p{{p'pf + 12 p'p p'^p + 2Q{p'2p?) \ ^ 
= 4Gf{pf + 2p2f) +1^ -9G'M.y ^ -^G'^f)' ) 

^ {4>'iF + G'p-2p'ip){8p'2p + 2p'p) ^ {p'p + 2p'2p){2p'2F-p'pf ^ 

— 9p2F ^P2F 

Since (p'^p + 2G'pp'^p < \G'^ {p'p - 2p'2pf , therefore 
X%{s*)- min A"/(s) 

> {(p'iF + G'F-2p'iF){^P'2F + '2p'F) ^ ^P'^ + ^-^2^^ ^'^^'^^ ~ ^^'^^^ 



^ \G'^{p'p-2p'^pf{Sp'2F + 2p'p) ^ ^^ Ap'F + 2p'2p){2p'2p-p'pf 



9p'2p 9P2F 
"J ^P2F 

Hence, when s* exists and the assumption (jS.Sp holds, we have 

min A^(s) = min{A^(0), A^(s*), A^(4)} = A^(4) < min A^^(s). 

0<s<4 0<s<4 

When s* does not exist, we can immediately get that 

min A^(s) = min{A^(0),A^(4)} = A^(4) < min Xf{s). 

We now combine this result with (j5.7p and summarize the stability relation between the fully 
atomistic model and the reconstruction-based local model by the following theorem. 
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Theorem 5.2. The relation between the stability of the reconstruction-based local model and the 
atomistic model depends on the sign of (f)2p + G'p ■ Ip'^p and is given by 

min A^(5) =A^(4) < min A^^(.) tf cj^'^p + G'p ■ 2p'^p > 0, 

0<s<4 0<s<4 

min A^(s)=mm{A^(0),A^(s*),A?.(4)}= min A^^^ ^/ cj^'^^ + G'p ■ 2p'^p = 0, (5.9) 

0<s<4 0<s<4 ' 

min XUs) = mm{A^(0), A^(5*), A^(4)} > min A^^ (s) ^/ c^'^p + G'p • 2p'^p < 0. 

0<s<4 0<s<4 

We note from the theorem that the reconstruction-based local model can be less stable than the 
fully atomistic model, which might cause stability problems when constructing a coupling method. 

6. Conclusion. 

In this paper, we give precise estimates for the lattice stability of atomistic chains modeled by the 
fully atomistic EAM model and the volume-based and the reconstruction-based local approxima- 
tions. We identify the critical assumptions for the pair potential, the electron density function, and 
the embedding function to study lattice stability. We find that both the volume-based local model 
and the reconstruction-based local model can give 0(1) errors for the critical strain. The critical 
strain predicted by the volume-based model is always larger than that predicted by the atomistic 
model, but the critical strain for reconstruction-based models can be either larger or smaller than 
that predicted by the atomistic model. 

Further research is needed to determine the significance of these results for multidimensional 
lattice stability and for atomistic-to-continuum coupling methods that couple an atomistic region 
with a volume-based local region through a reconstruction-based local region [5lll8j. 
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